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Abstract 

Dielectronic recombination (DR) of Na-like Fe 15+ forming Mg-like Fe 14+ via excitation of a 
21 core electron has been investigated. We find that configuration interaction (CI) between DR 
resonances with different captured electron principal quantum numbers n can lead to a significant 
reduction in resonance strengths for n > 5. Previous theoretical work for this system has not 
considered this form of CI. Including it accounts for most of the discrepancy between previous 
theoretical and experimental results. 

PACS numbers: 34.80. Lx 
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I. INTRODUCTION 



Understanding the properties of astrophysical and laboratory plasmas necessitates know- 
ing the ionization balance of the observed or modeled sources. This in turn depends on the 
underlying recombination and ionization processes. Of particularly importance are data for 
the electron-ion recombination process known as dielectronic recombination (DR) which is 
the dominant recombination mechanism for most ions in atomic plasmas jl[ 0] . 

The DR process can be expressed as 

e" + Af <— > [A^~ 1)+ ]** -> [A$T 1)+ ]* + u. (1) 

DR is a two-step recombination process which begins when a free electron e~ collides with 
an ion of element A with charge q+ and in initial state i. The incident electron collisionally 
excites a core electron of the ion with principal quantum number n c and is simultaneously 
captured forming a system of state j. This process is known as dielectronic capture. We use 
the word "core" here to distinguish initially bound electrons from the captured electron. The 
energy of the intermediate system [A^ ^ + ]** lies in the continuum and it may autoionize. 
DR occurs when the state j radiatively decays to a state / emitting a photon of angular 
frequency oo. This reduces the total energy of the recombined system to below its ionization 
threshold. Conservation of energy requires that the energy of the initial free electron and 
unrecombined ion balance that of the intermediate recombined system. Thus the relative 
kinetic energy of the incident electron equals the excitation energy AE of the core electron 
in the recombined system in the presence of a captured electron plus the binding energy Ej, 
of this captured electron in the recombined system, i.e., AE = E^ + E^. Because AE and 
Eb are quantized, E\. is quantized, making DR a resonance process. 

Here we explore a particularly nagging discrepancy between theory and experiment for 
the simple M-shell ion Na-like Fe 15+ forming Mg-like Fe 14+ . Good agreement between ex- 
periment {3] and theory has been found for Fe 15+ (ls 2 2s 2 2p 6 3s) DR via An c = and 1 
excitations of a 3s electron For DR via An r = 1 core excitation of a 21 electron, 

previous theoretical work has shown the importance of configuration interaction (CI) within 
a 2s 2 2p 5 3l3l'nl" complex for a fixed n [6:]. Including this single- n CI reduced the predicted 
resonance structure by a factor of 2. However, that work plus other recent work {4], s|, 
which also consider CI only within the same n complex, are still about a factor of 2 times 



larger than experiment [3!] for resonance energies above 650 eV [3] • These resonances involve 
captured electron quantum numbers of n > 5. 

In this work we investigate the cause of this discrepancy. We use the FAC (Flexible 
Atomic Code) 8j which is fully relativistic and utilizes the distorted wave approximation. 
We have made a more complete accounting of possible autoionization and radiative decay 
channels than previous theoretical works. Additionally, we pay particular attention to the 
effect of CI between different n configurations. This multi-n form of CI has been neglected 
in previous theoretical studies for this system. 

The rest of this paper is organized as follows. In Sec. [Ill we review the standard theoret- 
ical approach to calculate DR, discuss the autoionization and radiative decay channels we 
considered for Fe 15+ DR, and outline our approach to handling CI between different n com- 
plexes. We compare our theoretical results to experiment and previous theory in Sec. IHIl 
Lastly, we summarize our results in Sec. HVl 



II. THEORETICAL METHOD 



A. Standard approach 



We calculated DR using the independent process, isolated resonance (IPIR) approxi- 
mation [9|. This method treats radiative recombination and DR separately and neglects 
quantum mechanical interferences between the two and between DR resonances. These in- 
terference effects have been shown to be small in general [9j]. The DR cross section in the 
IPIR approximation for a multiply-excited intermediate state \<ftj) with resonance energy Ej 
is given to lowest order in perturbation theory by jlo| 



(^IDI^X^-IVI^) 2 



E — Ej + iTj/2 



(2) 



Atomic units are used here and throughout the paper unless otherwise noted. E is the 
collision energy and k is the linear momentum of incident free electron, both given in the 
electron-ion center-of-mass frame, g^ and gj are statistical weights, |^) is an initial recom- 
bining state which includes the incident free electron, and |$/) is a final bound state. The 
incident free electron is not correlated with the target ion. D is the dipole radiation field 



interaction 

N ' s=l 

where c is the light velocity, iV is the number of bound electrons before dielectronic capture 
and r s is the position vector of electron s from the nucleus. V is the electrostatic interaction 
between the N initially bound electrons and the continuum N + 1 electron 

N 1 



8=1 



r s - r N+ i 



The total resonance width 1^- is given by 

* / 

Aj k is the autoionization rate from j to any state k of A q+ and can be expressed as 

A% = 2n\(^ k \V\^)\ 2 . (6) 
A r jjr is the radiative decay rate from j to / which can be written as 

Aj / = K* / W i >| 2 . (7) 



The energy integrated cross section (i.e., resonance strength) of state j is given by 11] 



7i 2 g A A^Zf^jf , s 

in the approximation Tj <C Ej. The resonance strength can be re- written as the product of 
the dielectronic capture (DC) strength 

which is related to the autoionization rate through detailed balance, and the branching ratio 

Ba = - ^ £ - . (10) 

J a a i Ar \ ' 



B. DR channels of Fe 15+ 



For Fe 15+ DR via An c = 1 core excitation of a 21 electron, we considered the autoioniza- 
tion and radiative decay channels 
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e" + 2s 2 2p°3s *S 1/2 -> 2l 7 3l'3l"nV" ->■ i 




(11) 



(12) 



where Z < 1, Z' and /" < 2, and < 5. This includes the 2l 8 3l'3l" radiative decay channel 
which was not considered by [$| and the 2l 7 3l'3l" autoionization channel which was not 
included by [6j]. We also considered CI for all possible 2l 7 3l'3l" core configurations. Thus, 
unlike the previous theoretical work, 2s — > 31 promotions are included. 

For n > 6, the 2l 8 3l'nl"' configuration lies in the continuum and radiative decays to 
autoionizing levels are possible. These can then autoionize or radiatively stabilize via 



2l 8 3l'nl"' ->■ |2/ 8 3/"} +e 

2Z 8 3Z"n/'" 
2Z 8 3Z'3/" 
2l s 3l'n'l" 



->■ < 



> + oj. 



(13) 
(14) 



The branching ratio for these radiative decays to autoionizing levels followed by radiative 
cascades (DAC) can be given by 12 ] 



J2 t ^ jt + Zt>A r Jt ,B t , 



V 4" V 4' ' 

where the final states t and t' are below and above the ionization threshold, respectively. B t i 
is the branching ratio for radiative stabilization of t' and can be determined by evaluating 
Bj iteratively. 



C. Configuration interaction between different n resonance complexes 

We performed a large scale CI calculation between all 2l 7 3l'3l"nV" complexes from n = 3 
to 14. This allows us to consider CI between resonances with different captured electron 
principal quantum numbers. A large orbital sensitivity of DR to the choice of the initial 
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radial wave function has been reported in Mg 2+ 13]. For Fe 15+ this sensitivity is expected 
to be insignificant as a result of the high q of the ion. We explicitly investigated here the 
effects of optimizing radial wave functions on the 2l 7 3l'3l" and 2/ 8 3/' configurations of the 
recombining ion as well as on the 2l 7 3l'3l"3V" and 2l 8 3l'3l" configurations of the recombined 
ion. Only small differences in resonance strengths and energies were seen. In the end, radial 
wave functions were optimized on the 2l 8 3l'3l" configuration as that gave best agreement 
with the experimental results. 

The j CI mixed state 0™ for an n complex can be expanded in the f unmixed basis 0™,' 
of an n' complex using 



where c n iy denotes the mixing coefficient for the <f>^, basis. We calculated autoionization 
and radiative decay rates from the wave functions obtained using this CI mixing. Past 
studies have not considered CI mixing between different n complexes. In those studies 
autoionization channels of the form 2l 7 3l'3l"nV" — > 2l s nl"' + e~ and radiative decay channels 
of the form 2l 7 3l'3l"nl"' — > 2l s 3l'nl"' + uj were possible only between the states of same n. 
However, taking into account CI mixing between different n resonance complexes allows for 
additional autoionization and radiative decay channels. 



III. RESULTS 



A. Experiment 



Theoretical studies of Fe 15+ DR have been aided greatly by the merged-beams experi- 
mental results of shown in Fig. [TJ The measured data represent the DR cross section a 
times the relative collision velocity v T convolved with the experimental 



energy distribution yielding a rate coefficient (crv T ) The energy distribution is de- 

scribed by a flattened Maxwellian with a temperature of feeTy =2.4 meV along the beams 
and a temperature of ksT± = 0.1 eV perpendicular to the beams. Field ionization in the 
experiment limits the measured data to n < 86. 
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FIG. 1. DR resonance structure of Fe 15+ via An c = 1 core excitation of a 21 electron. The blue 
line shows the experimental results of |3j. The black line shows our results including CI only within 
the same n resonance complex. The red line shows our results including CI between different n 
resonance complexes for 3 < n < 14. See text for details. 



B. CI within the same n complex 



We performed explicit calculations of autoionization and radiative decay rates up to 
n = 14 and extrapolated for n from 15 to the experimental cutoff of 86. A simple hydrogenic 
scaling law was used for the resonance energies, the autoionization rates, and the radiative 
decay rates of the captured electron for n > 15. The radiative decay rate of the core electron 
was set to the n — 14 value for all n > 15. The calculated DR strengths were multiplied 
by v r and convolved with the experimental energy distribution of [3] . The results are shown 
in Fig. [TJ In the figure we have also labeled some of the strong resonances based on the 

7 



TABLE I. DR resonance energy and strength for the strongest level in each n complex for single-n 
CI. Also listed are the corresponding data for multi-n CI where 3 < n < 14. J denotes the total 
angular momentum of each level. For the level description, relativistically closed shells with J = 0, 
such as 2s 2 and 2po/ 2 , are omitted for brevity. 



Resonance energy (eV) Resonance Strength (10 19 cm 2 eV) 



Level 


J 


Single-n CI 


Multi-n CI 


Single-n CI 


Multi-n CI 


[(2p 1/2 3s)i3»3 /2 ]i/23d5/2 


3 


385.17 


383.80 


1.53 


1.54 


[(2p 1/2 3s)i3d3 / 2]i/24d 5 / 2 


3 


612.54 


612.51 


2.03 


1.94 


[(2pi/23s)i3d 3 / 2 ]i/ 2 5/7/2 


4 


694.94 


695.20 


1.64 


0.07 


lYo o \ o J 1 a J 

[(2pi/23s)i3d 3 /2ji/26d 5 /2 


3 


Ton C\ A 

730.04 


730.17 


1.33 


0.22 


[(2pi/23s)i3d 3/ 2]i/27d 5/ 2 


3 


754.21 


754.16 


0.96 


0.47 


[(2pi/23s)i3d 3/2 ]i/28rf 5/ 2 


3 


769.73 


769.75 


0.68 


0.40 


[(2pi/23s)i3d 3 / 2 ]i/ 2 94/2 


3 


780.30 


780.30 


0.53 


0.32 


[(2pi/ 2 3s)i3d 3 / 2 ]i/ 2 10d 5/2 


3 


787.81 


787.83 


0.41 


0.25 


[(2pi/23s)i3d3/ 2 ] 1 / 2 lld 5 /2 


3 


793.35 


793.35 


0.32 


0.21 


[(2pi/ 2 3s)i3d 3 / 2 ] 1/2 12d5/ 2 


3 


797.54 


797.55 


0.25 


0.13 


[(2pi/ 2 3s)i3d 3 / 2 ]i/ 2 13d 5 / 2 


3 


800.08 


800.81 


0.20 


0.13 


[(2pi/ 2 3s)i3d 3 / 2 ]i/ 2 14d 5/ 2 


3 


803.37 


803.41 


0.17 


0.07 



results of our calculations. Table Q] lists the calculated resonance energies and strengths 
for the strongest resonance level in each n complex for 3 < n < 14. Including the 2s — > 3/ 
promotion channel gives improved agreement between theory and experiment in the collision 
energy range of 400-500 eV. The resonances between 400-450 eV are in the better agreement 
with the experiment compared to the previous FAC results jlj] . Also the resonance at ~ 470 
eV does not appear unless this excitation channel is included. However including CI only 
within the same n complex does not remove the large discrepancy between theory and 
experiment for collision energies over 650 eV. 
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C. CI between different n complexes 



Explicit calculations for autoionization and radiative decay rates were again carried out to 
n = 14. For higher n, the extrapolation described in Sec. lIII Bl was performed. The convolved 
results are shown in Fig. [TJ Resonance strengths and energies are reported in Table H] for the 
selected resonances described in Sec. IIII B[ Figure [1] shows that above ~ 650 eV multi-n CI 
dramatically reduces the theoretical results compared to single-n CI. This reduction brings 
theory into very good agreement with experiment. The previous factor of 2 differences have 
been reduced to the level of tens of percent. The remaining differences near the series limit 



may be due to field ionization effects in the experiment as described by 14j . computational 
resources having limited the multi-n CI calculations ton < 14, or some combination thereof. 

A general sense for the importance of multi-n CI for n > 5 can be gained by looking at 
the mixing factors for the resonances listed in Table [B The mixing factor is given by 

Wn'l 2 = ^2\c n 'j'\ 2 , (17) 

r 

where the summation is over the f basis states in the n' complex. The mixing occurs between 
levels with the same parity, symmetry, and angular momentum. The mixing factors are 
plotted in Fig. [2j One sees that the n = 3 and 4 resonances of Table [I] are largely unmixed 
with other n' complexes but the n > 5 resonances can be strongly mixed. In particular the 
n = 5, 6, and 14 resonances are very strongly mixed with other n' complexes. 

To gain a more quantitative understanding on how multi-n CI can affect the predicted 
resonance strengths, it is helpful now to re-write Eq. (JS} using the expansion basis of Eq. (TIE)) 
which gives 

<j. — 71-2 9j Y^n'Y^jA C n'j^ 2 ^n'jHY^n'Y^j'(\^'j'\ 2 Y^f' / K'j'f) 
Ej J2n' Yjj> \ C n'j'\ 2 {J2k An'j'k + Yjf^n'j'f) 

Here A°:, 4 , 4 is the autoionization rate from the unmixed basis state to an initial state i 
and is given by 

A a nljH = 27r\my\<P]')\ 2 - (19) 

A®,j, k is given by Eq. (Tl9|) but changing i — > k. A r n ,-,j is the radiative decay rate from the 
to a state / and is given by 

Ab j ,f=\(*f\D\$)\ l . (20) 
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The coupling (i.e., interference) terms between different basis such as (^|V|0", ){(f>j„ IVI^jt) 
and flT)]^') {(J)™',', \T)\Q f) have been neglected just as in the 1P1R approximation. The 
dielectronic capture strength for the CI mixing can be re-expressed as 

5 - = ^|-EEi^'i 2 ^ ( 21 ) 

and the branching ratio for the CI mixing is given by 




3 4 5 6 7 8 9 10 11 12 13 14 



n' 

FIG. 2. Mixing factor \c n >\ 2 as a function of the mixing complex n' . Results are plotted for the 
resonances listed in Table [H Each curve is labeled by the initial n configuration before mixing 
between n' configurations is included. 

Now, taking the [(2pi/23s)i3<i3/2]i/26<i5/2 resonance level of the n = 6 complex listed in 
Table Has an example, we find that it mixes primarily with the basis levels listed in Table HTl 
Note that the autoionization rate A a nl -, i from the [(2p!/ 2 3s)i3(i3/2]i/26(i5/2 level to i is over a 
factor of 10 larger than the autionization rates from the other listed basis levels to i. This 
leads to a reduction in 5*dc by a factor of 6.5 when the values listed in Table [Til are used in 
Eq. (|2"Tj) . compared to what is calculated for the single-n CI case. 

The branching ratios for the listed levels which mix with the selected n = 6 resonance 
are all similar or smaller in value to that for this specific level. The resulting total branching 
ratio Bj given by Eq. (|22|) is reduced by a factor of 1.3 from the single-n CI value. Combining 
the multi-n values for Sbc an d Bj, we find a reduction for the resonance strength of the 

10 



TABLE II. Mixing basis level distribution for the [(2pi/23s)i3d 5 /2]i/2^d^/2 resonance listed in 
Table HI For the level description, relativistically closed shells with J = 0, such as 2s 2 , 2p 2 , 2 and 
2pg/ 2 > are omitted for brevity. The square of the mixing coefficient as defined in Eq. (|16|) is given 
in percent. A^,-^ is the autoionization rate from j' to i where i is the initial state 2s 2 2p 6 3s 2 5i/2 
of recombining ion, J2k^Ki'j'k + Ylf^-n'j'f * s ^ ne total autoionization and radiative decay rate of 
f, and Bj* is the branching ratio of f. Only basis levels where \c n 'j'\ 2 > 2% are listed. A total of 
10298 basis levels were included for this n = 6 resonance. 



j basis level 


1 |2 
l^j'n' \ 

(%) 


A a 

A j'n'i 
(s A ) 




+ E/ A j'n>f 

—Is 
) 


By 


[(2pi/ 2 3s)i3d 3 /2]i/26d 5 /2 


14.4 


2.64 x 10 13 


5.77 


x 10 13 


0.423 


[(2p3 /2 (J = 3/2)3s) 2 3p 3/2 ] 5/2 13/l 11/2 


8.4 


3.17 x 10 7 


5.24 


x 10 11 


0.032 


[(2p3 /2 (J = 3/2)3 S )23d 5/ 2] 1/2 6d 5/2 


6.8 


2.70 x 10 12 


8.46 


x 10 13 


0.051 


[(2p| /2 (J = 3/2)3a)i3p 3/2 ] 5/2 13^ 11/2 


5.8 


4.38 x 10 7 


8.24 


x 10 11 


0.012 


[(2pi/ 2 3f> 3 / 2 ) 2 3d 3 / 2 ] 3 / 2 5/ 5 / 2 


5.8 


1.40 x 10 10 


3.49 


x 10 13 


0.513 


(2p 3 /2 (J = 3/2)3p 2 /2 (J = 2)) 3/2 7r7 9/2 


4.8 


1.40 x 10 10 


1.43 


x 10 13 


0.003 


[(2s3p 1/2 )i3d 5/2 ]5/24s 


4.4 


4.22 x 10 10 


6.97 


x 10 12 


0.348 


[(2p3 /2 (J = 3/2)3a)i3d 5/2 ] 5/2 7s 


2.6 


4.22 x 10 10 


1.33 


x 10 13 


0.003 


[(2pi/ 2 3p 3 /2)l3d 3 /2] 3 /25/5/2 


2.3 


5.64 x 10 8 


2.08 


x 10 12 


0.131 


[(2s3p 1/2 )i3d 3/2 ]5/24s 


2.3 


1.00 x 10 11 


3.17 


x 10 13 


0.022 


[(2pi /2 3s)i3o? 3 /2]i/26(i5/2 level by a 


factor of £ 


.5, compared to the single-n CI results. 


This 



estimate agrees reasonably well with the factor of 6.5 reduction from the more complete 
calculation, as can be seen in Table [B The convolved DR resonance strengths for all n = 6 
resonances are displayed in Fig. |3j The reduction of the strong n = 6 resonances can 
be clearly seen. The reduction of the resonance strength for the other resonances listed 
in Table [I] can be similarly explained by multi-n CI mixing. In short, the DR resonance 
strength for strong resonances can be decreased by CI mixing with weak resonances. 

On the other hand, the DR resonance strengths for some weak resonances can be increased 
by CI mixing with strong resonances of different n' complexes. For example, the n = 13 
resonance strengths around 730 eV are largely enhanced by CI with the strong n = 6 
resonances at this energy, as is shown in Fig. |3j 
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FIG. 3. DR resonance structure of the n = 6 and n = 13 resonances for single-n CI and for multi-n 
CI. 

IV. SUMMARY 

We have demonstrated the importance of CI between resonances with different captured 
electron principal quantum numbers n for DR of Na-like Fe 15+ forming Mg-like Fe 14+ via 
An c = 1 core excitation of a 21 electron. Multi-n CI significantly reduces the theoretical 
resonance strengths for capture into n > 5 levels which overlap in energy with other many 
different n levels. This brings theory into very good agreement with experiment and removes 
a previously existing discrepancy between the two. The n = 4 levels are largely unaffected by 
multi-n CI because the energy separation between the n = 4 resonances and the interacting 
higher n resonances is large enough to render the multi-n CI unimportant. Such is not the 
case for the energy separation of the n > 5 resonances and those that they interact with, 
particularly for n = 5, 6, and 14. Additionally we have shown the importance of DR via 
2s — > 31 core promotions. 
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